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By means of QCD simulations on the lattice, we compute the coupling of the heavy-light mesons 
to a soft pion in the static heavy quark limit. The gauge field configurations used in this calculations 
include the effect of Nf — 2 dynamical Wilson quarks, while for the static quark propagator we 
use its improved form (so called HYP). On the basis of our results we obtain that the coupling 
g = 0.44 ± O.OSj^QQoi where the second error is fiat (not gaussian). 
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I. INTRODUCTION 



The static quark limit of QCD offers a simplified frame- 
work to solving the non-perturbative dynamics of light 
degrees of freedom in the heavy-light systems. That dy- 
namics is constrained by heavy quark symmetry (HQS): 
it is blind to the heavy quark flavor and its spin. As a 
result the total angular momentum of the light degrees 
of freedom becomes a good quantum number (jf ), and 
therefore the physical heavy-light mesons come in mass- 
degenerate doublets. In phenomenological applications 
the most interesting information involves the lowest lying 
doublet, the one with jf = (1/2)^, consisting of a pseu- 
doscalar and a vector meson, such as {Bq,B*) or {Dq,D*) 
states, where q € {m, d, s}. When studying any phe- 
nomenologically interesting quantity from the QCD sim- 
ulations on the lattice that includes heavy-light mesons 
(decay constants, various form factors, bag parameters 
and so on), one of the major sources of systematic uncer- 
tainty is related to the necessity to make chiral extrapola- 
tions. The reason is that the physical light quarks, which 
are expected to most significantly modify the structure 
of the QCD vacuum, are much lighter than the ones that 
are directly simulated on the lattice, rUg ^ rn^^d- Here 
by "q" we label the light quark masses that are attain- 
able from the lattice. Since the QCD dynamics with very 
light quarks is bound to be strongly affected by the ef- 
fects of spontaneous chiral symmetry breaking, a more 
suitable (theoretically more controllable) way to guide 
such extrapolations is by using the expressions derived 
in heavy meson chiral perturbation theory (HMChPT), 



which is an effective theory built on the combination 
of HQS and the spontaneous chiral symmetry breaking 
[SU{Nf)L ® SU{Nf)R -> SU{N})v]. Its Lagrangian is 
given by [l[ 

>Choavy = -tra Tr[i?aM)-Dft„Ht] + g tva Tt[H aH^,^ ^.A'^^-y^] , 



where 



Ha{v) 



(1) 



(2) 



is the heavy meson doublet field containing the pseu- 
doscalar, P"(w), and the vector meson field, P* "(w). In 
the above formulas, the indices a, h run over the light 
quark flavors, ^ = exp(i$//), with (f> being the matrix 
of {N'j — 1) pseudo-Goldstone bosons, and "/" is the pion 
decay constant in the chiral limit. We see that the term 
connecting the Goldstone boson (A^) with the heavy- 
meson doublet [i?(f )] is proportional to the coupling 
which will therefore enter into every expression related 
to physics of heavy-light mesons with = (l/2)~ when 
the chiral loop corrections are included. ^ Being the pa- 
rameter of effective theory, its value cannot be predicted 
but should be fixed in some other way. It can be re- 
lated to the measured decay width v\d* Dtt) 0], 
with the resulting value gcharm = 0.61(7). That value 
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^ A special attention should be given to the problem related to 
the presence of the nearby excited states as discussed in ref. 0|. 
Any precision lattice calculation cannot be fully trusted if the 
chiral extrapolations are made without discussing the problem 
of discerning the mixing with the = (1/2)+ states in the 
chiral loop diagrams. 
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turned out to be much larger than predicted by all of the 
QCD sum rule approaches but consistent with some 
model predictions such as the one in ref. in which 
a more detailed list of predictions with their references 
can be found. The large value for (7£)*£)7r-coupling was 
confirmed by the quenched lattice QCD study in ref. 
and recently also in the unqucnchcd case Since the 
charm quark is not very heavy, the use oi Qd'Ott to fix 
the value of g-coupling, via 



9 



-fir : 



(3) 



and its use in chiral extrapolations of the quantities 
relevant to i3-physics phenomenology may be danger- 
ous mainly because of the potentially large (^(l/m")- 
corrections. Unfortunately the decay B* Bit is kine- 
matically forbidden and therefore, to determine the size 
of g, we have to resort to a non-perturbative approach to 
QCD. Unlike for the computation of the heavy-to-light 
form factors, QCD sum rules proved to be inadequate 
when computing go'Din most likely because of the use 
of double dispersion relations when the radial excitations 
should be explicitly included in the analysis, as claimed in 
ref. Q. In this paper, instead, we compute the g-coupling 
on the lattice by using the unquenched gauge field con- 
figurations, with Nf = 2 dynamical light quarks, and 
in the static heavy quark limit. The attempts to com- 
pute this coupling in this limit were made in ref. Q , and 
very recently in ref. [lol |. On the basis of the currently 
available information, the coupling g in the static limit 
is indeed smaller than the one obtained in the charmed 
heavy quark case. 

In the remainder of this letter we will briefly describe 
the standard strategy to compute this coupling, list the 
correlation functions that are being computed to extract 
the bare coupling gq, as well as the axial vector renormal- 
ization constants. We then give details concerning the 
gauge field configurations used in this work, and present 
our results. 



II. DEFINITIONS AND CORRELATION 
FUNCTIONS TO BE COMPUTED 

In the limit in which the heavy quark is infinitely heavy 
and the light quarks massless, the axial coupling of the 
charged pion to the lowest lying doublet of heavy-light 
mesons, g, is defined via Q 



{B\A\B*{e))=gex. 



(4) 



where the non-relativistic normalisation of states \B^*^) is 
assumed, {Ba{v)\Bi,{v')) = Sab^iv — v'). For the heavy- 
light hadrons at rest (w = ■(/ = 0), the soft pion that 
couples to the axial current, = uj^'y^d, is at rest too, 
|(?| = 0. is the polarisation of the vector static- light 
meson. In the typical situation on the lattice we are 
away from the chiral limit {g — > gq), and the coupling 



(jq becomes the axial form factor whose value should be 
extrapolated to the chiral limit, in which the soft pion 
theorem relating the matrix element of the axial current 
to the pionic coupling applies @. 

The standard strategy to compute the above matrix 
element on the lattice consists in evaluating the following 
correlation functions: 



C2{t) 



C3{ty,tx) 



(ETr 



i^M/°755„,d(0,a;)75 



(5) 



1 + 70 



i,x,y 



W^jiSuiy, x)^a^Sd{x, 0)75 



where {■■■)„ denotes the average over independent gauge 
field configurations, the interpolating fields are P = h'y^q^ 
Vi = hjiq, with h{x) and q{x) the static heavy and the 
light quark field, respectively. In what follows, we drop 
the dependence on ty. In practice its value is fixed to 
one or several values as it will be specified in the text. 
In cq. ([5]) we also expressed the correlation functions in 
terms of quark propagators: the light ones, Sq{x, y), and 
the static heavy one, which becomes a Wilson line. 



Wy = 6ix-y) n 



x). 



(6) 



The latter is merely obtained from the discretized static 
heavy quark action [ll[ 



-Chqet = E ^^^^ 



6)t/i(x-6) ,(7) 



where for t/™^' ', the time component of the link variable, 
we use its improved form, obtained after applying the 
hyper-cubic blocking procedure on the original link vari- 
able, with the parameters optimized in a way described 
in ref. [H], namely with a = (0.75, 0.6, 0.3). That step 
is essential as it ensures the exponential improvement of 
the signal to noise ratio in the correlation functions with 
respect to what is obtained by using the simple product 
of link variables [l^ . 

The spectral decomposition of the three point function, 
given in eq. (O, reads 



where the sum includes not only the ground states (m = 
n = 0) but also their radial excitations (to, n > 0), which 
are heavier and thus exponentially suppressed. Note a 
shorthand notation, Z„ = \{0\h^^r,q\Bq)\, and the fact 
that we do not distinguish from couplings to the vec- 
tor interpolating operator because of the HQS. If the non- 
diagonal terms in the above sum were important (n ^ m) 
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FIG. 1: Ratio R{tx) in eq. i)J2)) as obtained from our data for 
all three sets and for which Kvai = «:soa ~ ^2 , with K2 specified 
in tabJeU This plot also shows the flatness of the signal of 
C3{tx) defined in eq. For completeness, we also note that 



the correlation function C^it^) would exhibit some expo- 
nential dependence in t^. In practice, it appears that the 
correlation functions C3{tx), as defined in eq. ([5]) are very 
flat (ta^-independent) for all the data sets that we use in 
this work and the details of which will be given in the 
next section (c.f. fig. This observation in fact agrees 
with what one can deduce from various quark models, 
and in particular from the one in ref. We will there- 
fore discard the non-diagonal terms in the spectral de- 
composition of C3{tx). We are still left with the problem 
of contamination of the desired signal (n = 0) by the 
axial transitions among radial excitations, n = m > 0. 
To solve that problem we should employ some smearing 
procedure and suppress the couplings of the source oper- 
ators to the radial excitations. To that purp ose we use 
the smearing technique proposed in ref. [14| . which es- 
sentially means that -in eq. ([5])- the interpolating fields 
are replaced by h{x)jc^q{x) h{x)j5q^{x), and similarly 
for the source of the heavy- light vector mesons, where 



r— k—x,y.z 



q(x^rk) X 

r 



,(8) 



and ip{r) = e^'^/^'{r + 1/2)^. The link variables on the 
right hand side of eq. ^ are fuzzed as discussed in ref. . 
After several trials we chose the smearing parameters to 
be i? = 1.3 and i?max = 4, to highly enhance the overlap 
with ground states. From the fits of our two-point func- 
tions computed with both the local ("loc") and smeared 



sources ("sm.") to two exponentials on the large inter- 
val 4 < t < 15, we obtain that Z™7Z^°'=- > 45, while 
^smy^ioc. ^ 0.05. More importantly, Zf^ /Zl"^- < 0.04, 
or it cannot be fitted, when it is completely absent. We 
therefore deduce that our smearing is efficient and the 
contribution of the radial excitations is most probably 
negligible. To further check this point we reorganized the 
operators in C3{tx) and fixed the transition operator [Ai 
in eq. (O] at t = 0, one source operator at ty = tfi^ = ~5, 
and have let the other source operator free (c.f. also 
ref. [III). In that situation the spectral decomposition 
looks as follows, 

C'.it^) ^Y.^le-'" (*-*-)(i?„|A,|i?:)£(") , (9) 



which allows us to check whether or not its effective bind- 
ing energy, with the smeared source operators. 



^Isitx) = log 



C'^itx + 1) 



(10) 



agrees with what is obtained from the two-point correla- 
tion functions. 



Cfr(0 = log 



C2{t) 
C2{t+l) 



(11) 



This is illustrated in fig.[H which we find satisfactory. Af- 
ter these checks, we extract ijq from the fit to a constant 
of the ratio, ^ 



1 



C3(i.) 



3 (zr-)^ 



(12) 



All our fits are made on the common interval, 5 < f^, < 8. 
On one ensemble of our gauge-field configurations we also 
checked that the value of (jq extracted from eq. is 
fully consistent with what is obtained if the computation 
is organized as in eq. 



A. Axial current renormalization constant 

The final ingredient necessary to relate the results of 
our calculation to the continuum limit is the appropri- 
ate axial current renormalization. We prefer to apply 
the same procedure to all our data sets and determine 
non-perturbatively the axial renormalization constant . 
To avoid any notational ambiguity we stress that in this 
subsection we discuss only the light bilinear quark non- 
singlet operators, P{x) = q{x)^5q{x), = q{x)^f^q{x), 
= q{x)jfj^j5q{x), i.e. no reference to the static 



^ The use of index "q" in gq should not be confusing to the reader. 
Here it simply labels the light quark directly accessed from our 
lattices. 
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2m(^^iP(z)K,(x)P(0)>,(14) 



FIG. 2: Comparison of the effective binding energy extracted 
from C2{t) and from C3(ta;) for the Set 2 and ki (c.f. table\^. 
Notice that the separation of the sources t = tx ~ tflx and in 
our case tfix = —5. 



Action 


P [a (fm)] # meas. ^sea ref. 


WP/ Wilson 
Set 1 


5.8 [0.054(2)] 50 ki = 0.1535 [IT] 
50 K2 = 0.1538 
50 K3 = 0.1540 
50 K4 = 0.1541 


Iwasaki/Clover 
Set 2 


2.1 [0.099(0)] 100 Ki = 0.1357 [18] 
100 K2 = 0.1367 
100 K3 = 0.1374 
100 K4 = 0.1382 


WP/Clover 
Set 3 


5.29 [0.075(1)] 60 ki = 0.1355 [19J 
80 K2 = 0.1359 
100 K3 = 0.1362 



TABLE I: Basic information on the sets of unquenched 
gauge field configurations with Nf — 2 dynamical Wilson 
quarks, the hopping parameters of which are specified for 
each set. "WP" stands for the Wilson-Plaquette gauge ac- 
tion, and Clover is a standard distinction to indicate the non- 
perturbative 0{a)-improved Wilson quark action. More in- 
formation on each set of configurations can be found in the 
quoted references. All lattice volumes are 24^ X 48. 



heavy quark will be needed in this subsection. To eval- 
uate ZAigi) ^se the hadronic Ward identity p^ . 
which is readily derived by imposing the invariance un- 
der the axial chiral rotations of {J2x^ii^)^ii^)) ^ ^^"^ 
{'ExMx)P{0)). One then obtains 



-Zy J d4z(^2,n|^4p(z)K,(x)A,(0)),(13) 

^ X 



where the integration 
not) include zero. 



volume V 
The bare 



(V) does 
quark 

2m 



(does 
mass de- 



(0) 



fined via the axial Ward identity reads, ^//i^^j 
{J2xAoix)P{0))/{J2sPix)PiO))- For notational sim- 
plicity in the above Ward identities we wrote Zy.A = 
Zv,A{gl,amq). 



III. LATTICE DETAILS AND RESULTS 

We use the publicly available gauge field configurations 
generated with Nf ~ 2 dynamical light ("sea") quarks 
which were produced by using the Wilson gauge and Wil- 
son quark actions. In tablc[T]we provide a basic informa- 
tion on the data sets used in this letter. Concerning the 
discretizcd Yang-Mills part, the configurations explored 
in this letter were generated by the standard Wilson pla- 
quette action and (in one of the sets) by its improved 
form, known as the Iwasaki action. The effects of dy- 
namical quarks in the QCD vacuum fiuctuations are sim- 
ulated by using the Wilson quark action, both the ordi- 
nary one, and its non-perturbatively C'(a)-improved ver- 
sion, which is usually referred to as the "Clover" -action. 
From the publicly available configurations we chose those 
with small lattice spacings, a < 0.1 fm. In table[I]we also 
provide the references containing detailed information 
about the simulation parameters and the algorithms used 
in producing these configurations. The values of lattice 
spacings, given m table m are obtained from rg/a, com- 
puted on each of these lattices, extrapolated to the chiral 
limit and then by choosing tq — 0.467 fm. Other popu- 
lar choice is rg = 0.5 fm, which would make the lattice 
spacing 7% larger. To our purpose that error on fixing 
the lattice spacing is, however, completely immaterial. 
We should emphasize that we do not work in the par- 
tially quenched situations. Instead, we fix the hopping 
parameter (Kg) of our valence light quark in correlation 
functions ([5]) and in those appearing in eqs. (|13ll4p to 
be equal to that of the corresponding dynamical ("sea") 
quark, also fisted in table HI 

In this paper we do not use the so-called "all-to-all" 
propagators. The feasibility study of using that tech- 
nique in the computation of ^-coupling has been made 
recently in ref. [lO[ , showing the substantial reduction in 
statistical errors. We plan to adopt that technique in our 
future studies. In table|TT]we provide the list of all results 
relevant to the subject of this letter, that we directly ex- 
tracted from the correlation functions computed on all 
lattices from table HI For an easier comparison, the val- 
ues of the pseudoscalar light meson masses, as wefi as of 
the bare light quark masses inferred from the axial vector 
Ward identity, are given in lattice units. They are fully 
consistent with those reported in refs. [13, [H, [l^. Con- 
cerning the renormalization constants Zy.A, they are ob- 
tained from the "light-light" correlation functions which 
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FIG. 3: Qq computed from the ratio in eq. i)J2D for all of 
our lattice data sets listed in table U after accounting for 
the axial current renormalization constants computed on the 
same ensembles of gauge field configurations. They are plot- 
ted as a function of the light pseudoscalar meson ("pion") 
mass squared (in GeV^j. 



we computed on the lattice and then combined to verify 
the Ward identities in eqs. (|13ll4p . After inspection, we 
found the common plateau-region for all our 11 data sets 
to be between 10 < t < 14. Finally, in the three-point 
correlation function in cq. ([S]) the fixed source operator 
is set at = 13, and -as already mentioned- the results 
for gq are obtained from the lit to a constant in eq. (jl2p . 
on the interval 5 < < 8. We checked that our results 
remain stable when ty ~ 12. Directly extracted values 
for the bare couplings cjq, from all of the lattice data sets 
considered in this letter, i.e. before multiplying them by 
its corresponding Za, are listed in table HIl In fig. El in- 
stead, we plot the renormalized coupling gq, as a function 
of the squared light-light pseudoscalar meson ( "pion" ) , 
mass now given in physical units (in GeV'^). That conver- 
sion is made by computing tqitIt^ for each of our data sets 
and then use rp = 0.467 fm (or, = 2.367 GeV"^). We 
reiterate that opting for tq = 0.5 fm [tq = 2.534 GeV~^), 
does not alter our final results in any significant way. 

The last step to reach the coupling which is our final 
goal, is to make the extrapolation to the chiral limit. To 
that end we attempt either a simple linear fit or a fit 
guided by the expression derived in HMChPT [2^, i.e.. 



5^ =50 



1 



(47r/)- 



:ml log(m^) + coml 



(15) 
(16) 



where g^ is then the soft pion coupling that is to be used 
in applying the HMChPT formulas when extrapolating 
the phcnomcnologically interesting quantities computed 
on the lattice to the physical light quark mass limit. From 



fig. [3] it is obvious that this task is quite difficult if one is 
doing it separately for each (3. More specifically, applying 
the linear fit (fTSll to each of our data sets we obtain 



= {(0.40 ± 0.15)i, (0.52 ± 0.07)2, (0.54 ± 0.06)31(17) 

while from the fit to HMChPT (HH) we get 

go = {(0.36 ± 0.11)i, (0.43 ± 0.04)2, (0.48 ± 0.04)3} j(18) 

where the index on the right hand side labels the data 
sets like in the tables U and [Til As it could have been 
anticipated from eq. (fTHl) . the results of the HMChPT 
fit (50) are lower than the results of linear extrapolation 
(5iin)- The values obtained from different sets are consis- 
tent within the errors. It is obvious that we cannot make 
a precision determination of this coupling yet, but it is 
clear that the unquenched lattice data also point to the 
fact that the g coupling is considerably smaller in the 
static-heavy quark limit than in the case of the heavy 
charm quark. Since the heavy quark is only a spectator, 
this information -that the 1 / m^-corrcctions arc large- is 
significant, and somewhat surprising. If one simply feeds 
the difference by a linear 1 / mc-tcrm, it is quite interest- 
ing to notice that from the light cone QCD sum rules 
one get a similar size is such a correction in spite of the 
fact that the absolute value for the pionic couplings were 
considerably underestimated [1] . 

Since we are not aiming at a percent-level precision 
determination of this coupling we can try and see what 
happens if all the data are combined and fit them to- 
gether to eqs. (jl5ll6p . We are, of course, aware that our 
three sets suffer from different discretization errors but 
since the lattice spacing is small (a < 0.1 fm) and the 
common renormalization procedure has been applied to 
all of them, it is reasonable to assume that the remaining 
discretization errors are not likely to matter, in view of 
our statistical error (~ 10%). If we combine all of our 
data, we then obtain 

giin = 0.51 ± 0.04 , Clin = (0.21 ± 0.12) GcV"^ , (19) 

while with the HMChPT formula ^ we have 

.90 = 0.44 ± 0.03 , Co = (0.40 ± 0.12) GeV"^ 



(20) 



Another possibility is to exclude the data with > 
0.6 GeV^ which gives 50 = 0.46 ± 0.04. We also checked 
that our resulting ^0 is insensitive to the variation of 
/ e (120, 132) MeV, latter being /p^J"^-. 

Before concluding we should compare our result to the 
existing unquenched value for g reported in ref . . The 
main advantage of the calculation presented in ref. 
with respect to ours is that they used the so-called all- 
to-all light quark propagators so that their resulting sta- 
tistical errors are much smaller. However the lattices 
we used here are finer and the associated discretization 
errors should be smaller. In addition, here we also use 
various gauge and quark actions, to show that our re- 
sults are robust in that respect too (of course within our 
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error bars). A reasonable comparison with ref. [10[ can 
be made by using our results from Set-2 because these 
data correspond to the same gauge and quark actions 
as those used in ref. [l^l, although the lattice spacing 
we use here is smaller. Comparing the bare quanti- 
ties, we see that -for example- when the pion mass is 
0.75 GcV, from ref. [3 we read 5^=1-8 = 0.68(1), 
^^=1.95 ^ 0.69(1), while our g^=2-i = 0.69(6). There- 
fore they fully agree although our statistical errors are 
much larger. Using the perturbativc (boosted) 1-loop ex- 
pression (bpt), the result for the overall renormalization 
constant in all three cases are equal among themselves 
within less than 1%, so that the renormalization con- 
stants computed in bpt would not spoil this comparison, 
which seems to indicate that the discretization errors are 
indeed small. The step in which we go beyond ref. 
is that we evaluate the axial renormalization constant 
non-pcrturbatively, Z'^'^ . This is particularly important 
when using the data obtained with Iwasaki gauge action 
because in that case the strong coupling is very large and 
the use of perturbation theory is far from being justi- 
fied. Various boosting procedures can lead to various es- 
timates of Za- We show that the boosting procedure used 
in refs. [13, E3| leads to the values very close to our non- 
perturbative estimate. More precisely, at our /3 ~ 2.1, we 
have Z^VZ^P' = 0.90(1), 0.94(2), 0.97(3), 0.97(6), when 
going from the heaviest to the lightest quark mass. 

IV. CONCLUSIONS 

In this letter we report on the results of our calcula- 
tions of the soft pion coupling to the lowest lying doublet 
of static heavy-light mesons. From our computations, 
in which we use the fully unqucnchcd set-up and three 
different sets of gauge field configurations, all produced 
with Wilson gauge and fermion actions, we obtain that 
^0 = 0.44±0.03^o go- The second error reflects the uncer- 
tainty due to chiral extrapolation and it is the difference 
between the results of linear fit and the fit in which HM- 
ChPT is used. If our result is to be used in the chiral 
extrapolations of the phenomenologically relevant quan- 



tities in i?-physics the second error should be considered 
as flat. The reason is that our central value is obtained 
via HMChPT fit, but since the domain of applicability of 
HMChPT is still unclear Q -as of now- both results (ex- 
trapolated linearly or by using the chiral loop correction) 
are equally valid. 

On the more qualitative level, our results show/confirm 
that this coupling is smaller in the static limit than what 
is obtained when the heavy quark is propagating and is 
of the mass equal to that of the physical charm quark, 
g^i,^^^ = 0.68±0.07 0]. It is intriguing that the 0{l/m'^) 
corrections are quite large for the quantity in which the 
heavy quark contributes only as a spectator. That fea- 
ture can be safely studied on the lattice by means of the 
rclativistic heavy quark action of ref. [2l| , and we plan to 
do such a study. An obvious perspective concerning the 
determination of ^0 is to further reduce the errors, both 
statistical (by using the "all-to-all" propagator technique, 
like in ref. [10|), and the systematic ones (in particular 
those associated with chiral extrapolations) . Once a per- 
cent accuracy is reached, it will be important to study 
carefully the effects of mixing with the lowest heavy-light 
excited states [those with j[ = (1/2)+]. The expression 
derived in HMChPT which accounts for those effects in 
gq already exist (first paper in ref. [l^l), but their use re- 
quires the knowledge of another pionic coupling, the one 
that parametrizes the 5-wave pion emitted in the transi- 
tion from a (1/2)+ (1/2)^ states. Finally, the numer- 
ical tests concerning the impact of inclusion of heavier 
quarks in the vacuum fluctuations (s and c) on the size 
of go J would be highly welcome too. 



Acknowledgments 

We thank the SPQcdR, CP-PACS and QCDSF col- 
laborations for making their gauge field configurations 
publicly available, The supports of 'Flavianet' (EU con- 
tract MTRN-CT-2006-035482) and of the ANR (con- 
tract DIAM ANR-07-JCJC-0031) arc kindly acknowl- 
edged too. Laboratoire de Physique Thcorique is unite 
mixte de Recherche du CNRS - UMR 8627. 



[1] R. Casalbuoni, A. Deandrea, N. Di Bartolomeo, 

R. Gatto, F. Feruglio and G. NarduUi, Phys. Rept. 281 

(1997) 145 hep-ph/9605342^ . 
[2] D. Becirevi c, S. Fajfer and J. F. K amenik, JHEP 0706 

(2007) 003 [arXiv:hep-ph/0612224 ; 
[3] A. Anastassov et al. [CLEO Collaboration], Phys. Rev. 

D 65 (2002) 032003 hep-ex/0108043 . 
[4] A. Khodja mirian et al, Phys. Lett. B 457 (1999) 245 

[hep-ph/99 03421 , see also P. Ball and R. Zwicky Phys. 

Rev. D 71 (2005) 014015 hep-ph/0406232 . 
[5] D. Becirev ic and A. L. Yaouanc, JHEP 9903 (1999) 021 

'arXiv:hep-ph/9901431 . 
[6] A. Abada et al, Phys. Rev. D 66 (2002) 074504 



[hep-ph/020 6237] . 

[7] D. Becirevic and B. Haas. [arXw:090 3. 2407 [hep-lat]. 
[8] D. Becirevic, J. Charles, A. LeYaouanc, L. Oliver, 
O. Pene and J. C. R aynal, JHEP 0301 (2003) 009 
'arXiv:hep-ph/0 212T77] . 
[9] G. M. de Divitiis et al. [UKQCD Collaboration], JHEP 
9810 (1998) 010 hep-lat /9807032 ; A. Abada et al., 
JHEP 0402 (2004) 016 hep-lat/0310050 . 
[10] H. Ohki, H. Matsufuru and T. Onogi, Phys. Rev. D 77 

(2008) 094509 ^arXiv :0802.1563 [hep-lat]]. 
[11] E. Eichten and B. Hill, Phys. Lett. B 234 (1990) 511. 
[12] A. Hasenfratz and F. Knechth, Phys. Rev. D 64 (2001) 
034504 ,arXiv:hep-lat /0103029. . 



7 



[13] M. Delia Morte, S. Durr, J. Heitger, H. Molke, J. Rolf, 
A. Shindler and R. Sommer [ALPHA Collaboration], 
Phys. Lett. B 581 (2004) 93 [Erratum-ibid. B 612 (2005) [19] 
313] 'arXiv:hep -lat/03070 21 . 

P. Boyle [UKQC D"boIia boration] , J. Comput. Phys. 179 
(2002) 349 hep-lat/9903033 . [20] 
G. Martinelli, S. Petrarca, C. T. Sachrajda and 
A. Vladikas, Phys. Lett. B 311 (1993) 241 [Erratum- 
ibid. B 317 (1993) 660]; M. CrisafuUi, V. Lubicz 
and A. Vladikas, Eu r. Phys. J. C 4 (1998) 145 
arXiv:hep-lat /9707025] . T. Bhattacharya, R. Gupta, [21] 
W. J. Lee and S. R. Sharpe, Phys. Rev. D 63 (2001) 



[14] 
[15] 



and A. Le Yaouanc, 



074505 [arXiv:he p-Iat /0009038 
[16] D. Becirevic, E. Chang 

larXiv:0905.3"352l [hep-lat]. 
[17] D. Becirevic et al., N ucl. Phys. B 734 (2006) 138 

[arXiv:hep-lat/0510 0"T4] . 
[18] A. Ah Khan et al. [CP-PACS Collaboration], Phys. Rev. 



D 65 (2002) 054505 [Erratum-ibid. D 67 (2003) 059901] 
arXiv;hep-lat/0105015^. 

(2006) 

M 



A. A. Khan et al, P hys. 
"arXiv:hep-lat /0603028 



Rev. D 74 
Gockeler et al.. 



094508 
Phys. Rev. 

D 73 (2006) 054508 a r'Xiv:hep-lat/0601004 | . 
S. Fajfer and J. F. Kamenik, Phys. Rev. D 74 
(2006) 074023 arXiv:hep-ph/0606278 ; D. Becirevic, 

(2003) 
Nucl. 



S. Prelo vsek and J. Zupan, P hys 
054010 'arXiv:hep-lat/0210048 



Rev. D 67 
W. Stewart, 



Phys. B 529, 62 (1998)7 
A. X. El-Khadra, A. S. Kronfeld and P. B. Mackenzie, 
Phys. Rev. D 55 (1997) 3933 arXiv:hep-lat/9604004l; 
N. H. Christ, M. Li and H. W. Lin, Phys. Rev. D 
76 (2007) 074505 ,arXiv:hep-lat/0 608006 ; Y. Kayaba 
et al. [CP-PACS CoUab orationj, JHEP 0702 (2007) 019 
"arXiv:hep-lat /061 1033] . 



8 



p 






(0) 


Zv{go,mq) 




9q 


5.8 

Set 1 


0.1535 
0.1538 
0.1540 
0.1541 


0.262(4) 
0.236(4) 
0.221(3) 
0.182(7) 


0.0333(6) 
0.0260(3) 
0.0215(4) 
0.0180(4) 


0.720(38) 
0.705(41) 
0.643(57) 
0.654(61) 


0.908(50) 
0.858(51) 
0.795(70) 
0.827(81) 


0.683(52) 
0.714(71) 
0.742(81) 
0.628(62) 


2.1 

Set 2 


0.1357 
0.1367 
0.1374 
0.1382 


0.631(2) 
0.519(2) 
0.422(2) 
0.298(2) 


0.1078(5) 
0.0743(3) 
0.0513(4) 
0.0252(4) 


0.742(13) 
0.751(25) 
0.754(38) 
0.752(82) 


0.822(13) 
0.839(18) 
0.847(26) 
0.829(50) 


0.738(42) 
0.736(56) 
0.691(55) 
0.622(61) 


5.29 
Set 3 


0.1355 
0.1359 
0.1362 


0.327(2) 
0.245(2) 
0.155(2) 


0.0357(3) 
0.0206(2) 
0.0086(3) 


0.772(2) 
0.770(4) 
0.760(10) 


0.793(5) 
0.786(9) 
0.779(14) 


0.683(63) 
0.658(61) 
0.711(62) 



TABLE II: Direct numerical results extracted from the correlation functions calculated on all of the ensembles of the lattices with 
parameters enumerated in table U 



